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Exploiting previous results on Markovian dynamics and fluctuation theorems, we study the conse- 
quences of memory effects on single realizations of nonequilibrium processes within an open system 
approach. The entropy production along single trajectories for forward and backward processes is 
obtained with the help of a recently proposed classical-like non-Markovian stochastic unravelling, 
which is demonstrated to lead to a correction of the standard entropic fluctuation theorem. This 
correction is interpreted as resulting from the interplay between the information extracted from 
the system through measurements and the flow of information from the environment to the open 
system: Due to memory effects single realizations of a dynamical process are no longer independent, 
and their correlations fundamentally affect the behavior of entropy fluctuations. 
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I. INTRODUCTION 

Since many years fluctuations in far-from-equilibrium 
processes have been studied both in a classical [IH3] and 
in a quantum context j3HS]. The large interest they at- 
tract is due, for instance, to their deep connection with 
equilibrium thermodynamic quantities characterizing a 
system interacting with an environment [5], as well 
as with quantum information theory 9J. In the quan- 
tum realm the consequences of such an interaction on 
the fluctuations of physical quantities of the open sys- 
tem has been extensively studied only under the Marko- 
vian approximation and/or in the weak coupling limit 
[TTH ITT] . Different methods have been developed, allow- 
ing for the formulation of some theorems analogous to the 
well-known classical fluctuation theorems [T^lU4) . How- 
ever, an extension of these approaches to non-Markovian 
processes has only partially been attempted [15] and is far 
from being a settled problem. Taking into account non- 
Markovian dynamics in nonequilibrium thermodynamics 
is, however, an issue of great interest, since recently mem- 
ory effects have been extensively investigated, both from 
a theoretical [16H18) and from and experimental point of 
view |19) . On a more fundamental level, the possibility 
of characterizing a non-Markovian extension of thermo- 
dynamic concepts is very closely related to the investiga- 
tion on the quantum limits of thermodynamics [201 121j , 
recently attracting a great deal of attention. 

In this work we extend the method developed in 
Ref. |10j to the non-Markovian regime by employing 
recently constructed stochastic formulations of classical 
and quantum non-Markovian master equations [22H24j . 
In this way we are able to account for memory effects 
in far-from-equilibrium entropic fluctuations and thus to 
formulate a non-Markovian generalization of an entropic 
fluctuation theorem, highlighting also how and why it 
differs from its Markovian counterpart. 

The paper is structured as follows. In Sec.[TJjwe define 
our physical scenario and introduce its general dynamical 
description. SectionllTTlis a brief review of the Markovian 



approach of [TO], slightly modified to fit with our goals, 
and explicitly emphasizes the role of the structure of the 
associated master equation. In this way the reason for 
which such a method cannot account for non-Markovian 
processes will become clear. Section IV is devoted to 



the introduction of the main ideas and results of this 
work, namely a generalization of |10j which exploits a 
stochastic unravelling of non-Markovian processes. In 
Sec. [V] we comment on the physical interpretation of the 
non-Markovian fluctuations described by our approach, 
and present some final remarks and conclusions. 



II. THE MASTER EQUATION 

We consider an open quantum system with free Hamil- 
tonian i?s, interacting with an environment (free Hamil- 
tonian He) via an interaction term Hi. Let us suppose 
the free Hamiltonian of the open system to show a cer- 
tain time dependence caused by some parameters being 
externally modified in time. Such an external modifica- 
tion drives the system far from its initial equilibrium (or 
stationary) state and towards a final state characterized 
by different values of these control variables. During its 
evolution the system interacts with its environment, our 
goal being to study the effect of such an interaction on 
fluctuations of thermodynamic quantities. We employ 
a time-convolutionless (TCL) master equation [25] for 
the open system density matrix p{t), describing a gen- 
eral quantum dynamics: 

d -o{t)=K,{t)p{t), (1) 



df 



where 



K{t)p = -i[H(t),p] 

+ £ 7l (t) L (t) P A\ (t) - \ \A\ (t)A t (t) , p) 



(2) 



The Hamiltonian H(t) in Eq. ^ describes the unitary 
part of the open system evolution, which is given by its 
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free Hamiltonian Hs(t) plus a renormalization term due 
to the interaction with the bath. The nonunitary part 
of the evolution, describing dissipation and dephasing, is 
accounted for by a set of generally time-dependent Lind- 
blad operators Ajft) and corresponding relaxation rates 
7i(i). Equation (JTJ) describes Markovian as well as non- 
Markovian processes in terms of a master equation which 
is local in time |26) . 

The starting point of our analysis on nonequilibrium 
entropic fluctuations is the expression for thermodynamic 
ensemble quantities, and in particular for the time vari- 
ation of the von Neumann entropy along the nonequilib- 
rium dynamics itself. To this end, we write the state of 
our open quantum system as [TU] 



pit) = £>„(t)|a(f))(a(t)|, 



(3) 



where the set *B (t) = is a time-dependent or- 

thonormal basis of the open system Hilbert space instan- 
taneously diagonalizing p(t), and ^2 a p a it) = 1- Using 
the generator of the dynamics given in ([2| and calculat- 
ing the mean value of Eq. ^ for a state \b(t)) € Q5(t), 
one easily obtains a Pauli-type master equation for the 
evolution of the populations fJ, a (t), 



&(*) = J2 ( R ba{t)Ha{t) - Rab{t)Hb{t) 



(4) 



where the total instantaneous transition rate Rb a (t) be- 
tween two states \a(t)) and \b(t)) belonging to 23(£) is 
defined as 



R ba (t)=J27i(t)\{m\Mt)\a(t))f 



(5) 



These transition rates will turn out to be crucial in the 
expression for all thermodynamic quantities of interest. 
Moreover, it is their time behavior which we employ to 
characterize the occurrence of memory effects during the 
process: As will be discussed in more detail in Sec. IV 
the process described by the Pauli master equation (|4j) 
is defined to be non-Markovian if and only if at least 
one of the transition rates Rb a (t) between two particu- 
lar instantaneous eigenstates of p(t) temporarily becomes 
negative. 



A. Entropies 

Evaluating the von Neumann entropy in the instanta- 
neous eigenbasis of pit), we obtain its time derivative S 
in the form 



£(i) = -5^£i(t) ln Mb (£). 

b 

Using Eq. Q in Eq. ^ we find 



S(t) = -J2n a (t)R ba it) In 



^bjt) 



(6) 



(7) 



A similar expression was derived in Ref . |10j , but the use 
of a TCL master equation in our approach allows us to 
express transition rates explicitly in terms of Lindblad 
operators. This clarifies the physical framework we are 
working in and, as will become evident in Sec. |IV[ explic- 
itly shows where and how memory effects come into play 
in the case of non-Markovian dynamics. 

Having at our disposal the expression for the time 
derivative of entropy, and following usual prescriptions 
of nonequilibrium thermodynamics [27, 28], we can write 
Eq. ([7]) as a sum of two different contributions as Sit) = 
S e (t) + Si(t), having defined 



&(t) = -5>-(t)-Rfca(*)ta E! " :(,) 

a j b 

Siit) = Y,Va(t)Rba(t)ln 



Rabit) ' 
Vait)Rbait) 



a . b 



fJ.b(t)R a b(t) 



(8) 

(9) 



as, respectively, the entropy flux between system and en- 
vironment and the total entropy production. Equations 
([7]), Q and |9| describe the time dependence of entropy 
due to the ensemble dynamics described by a TCL mas- 
ter equation. The irreversibility of the process is char- 
acterized by a nonzero rate of entropy production Si(t) 
inside the system which, furthermore, in the case of a 
Markovian dynamics never becomes negative. It is worth 
stressing that the definition ^ for the entropy produc- 
tion coincides with the negative time derivative of the 
relative entropy of pit) and the stationary state p s t&t of 
the dynamics, provided the latter exists and detailed bal- 
ance holds. 



III. THE MARKOVIAN CASE 



B. Fluctuations 



In the case of a purely (and possibly time-dependent) 
Markovian dynamics, the rates given by Eq. ([5| never 
become negative. In what follows we will use indices a 
and b to label vectors in the instantaneous eigenbasis of 
pit) (recall Eq. ([3])), and an index i to label the possible 
decay channels described by the set of Lindblad opera- 
tors Aiit) in |2]). We will furthermore, for the sake of 
brevity, sometimes suppress the time arguments of the 
instantaneous eigenvectors and eigenvalues of pit). 



The above defined quantities characterize the physics 
of the open quantum system on an ensemble level. Such 
a picture, while allowing us to derive suitable expressions 
for many quantities of interest, lacks however a character- 
ization of single nonequilibrium processes and, in partic- 
ular, of their intrinsically fluctuating physical quantities. 
The goal of this section is to obtain a fluctuation theorem 
for the entropy production along single realizations of the 
ensemble dynamics. By definition, a fluctuation theorem 
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for a quantity Q characterizing a thermodynamic system 
is an expression for the ratio of the probability of such a 
quantity having the value q along a particular nonequi- 
librium process, and the probability of the same quantity 
having a value — q along the backward realization of the 
same process. In order to describe these fluctuations we 
employ the master equation Q which in the present case 
only involves positive transition rates and can thus be re- 
garded as a differential Chapman-Kolmogorov equation 
for a classical, Markovian stochastic jump process with, 
in general, time-dependent rates Rb a {t). Let us consider 
a particular, yet generic realization of this process, 



\a (t )) -> K(*i)) 



\a N (t N )), (10) 



consisting of N jumps at times ti,t2, ■ • • , tjy between well 
defined states belonging to the instantaneous eigenbases 
of p(t) . According to the master equation Q the proba- 
bility of the trajectory ( 10 ) can be written as 



3=0 
N-l 

x H R aj+iai (t j+1 )dt j+1 , (11) 
3=0 

where the first factor gives the probability for the system 
to start its trajectory in the state |ao), the second factor 
gives the probability that there are no jumps between 
the times tj and ij+i, and the third factor represents 
the probability of having N jumps within infinitesimal 
time intervals dtj around tj between the states of the 
trajectory (10 1. 



It is important to emphasize that the stochastic de- 
scription given by Eq. (11) is based on the Pauli-type 
master equation Q and thus correspond to the stan- 
dard stochastic unraveling of a classical Markovian mas- 
ter equation |29| , in which the probability for a transition 
from state \a(t)) to state \b(t)) during the time interval dt 
is determined by Rb a (t)dt with the rate given by Eq. ([5]). 
Thus we follow here the interpretation proposed in |10j 
to identify the fluctuations of single nonequilibrium pro- 
cesses with those described by the evolution equation 
Q for the populations of the density matrix. This in- 
terpretation and the underlying physical picture has to 
be carefully distinguished from the interpretation of the 
stochastic wave function methods (see, e.g., Ref. [30] and 
references therein) for the open system dynamics given 
by a quantum master equation of the form of Eq. in 
terms of a continuous measurement of the environment. 

The backward process corresponding to ( |10| is de- 
scribed by the trajectory 



|ajv(iiv)) -> |o/v-i(*/v-i)) -> 



M*o))- (12) 



The probability pb for such a backward process is defined 
analogously to what has been done for the forward one, 
taking into account the jumps in the sequence given in 
Eq. ( 12 ) such that the conditioned non-jump evolution 



probability is then the same for forward and for backward 
processes, while the jump rates are reversed, 



N-l 



^)n^ +ldrWT) 



Pb = 



3=0 
N-l 

x U R a . aj+1 (t j+1 )dt j+l , (13) 

Thus, we find that the logarithm of the ratio of the two 
probabilities takes the form 



ln P/ =ln ^o(*0) 



Pb > ON (*iv) 

The first term 



^ln |^±m±4 . (1.J) 



3=0 



R aj a J + 1 {tj + l) 



As = In 



Ha {to) 



(15) 



^a N {tN) 

represents the change of the von Neumann entropy, while 
N-l 



As e 



3=0 



Raj + iaj (tj + l) 
R aj a j+ i{tj+l) 



(16) 



yields the entropy flux integrated along a single trajec- 
tory. The average over all possible trajectories leads to 
the expressions ([7| and §8h, respectively. Defining 



a = As — As e 



(17) 



as the total entropy production along a single trajec- 
tory, we thus obtain an entropic fluctuation theorem for 
Markovian dynamics, 



Pb(-v) 



(18) 



from which the quantum analog of Crooks theorem |31j 
and the quantum Jarzynski equality 32 directly follow. 



IV. THE NON-MARKOVIAN CASE 

What happens if we introduce memory effects into the 
dynamics? The generator ([2| of the TCL master equa- 
tion describing a non-Markovian time evolution keeps the 
same structure as before, but the decay rates 7i(i) can 
temporarily become negative. This allows us to describe 
in full generality any non-Markovian quantum process, 
the TCL formulation being very general, requiring only 
the map to be invertible and differentiable with respect to 
time [33] . We are going to study the thermodynamic con- 
sequences of this behavior. An attempt to take into ac- 
count negative decay rates in the formulation of a fluctu- 
ation theorem has already been performed in 15 where, 
however, a clear formulation of a non-Markovian fluctua- 
tion theorem was not given. In particular the approach of 
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[T5] works well as long as the transition rates Rb a {t) stay 
positive, which however is the characterization we gave 
of a Markovian process. A study of nonequilibrium fluc- 
tuations, then, has not been performed yet for the class 
of processes which we define as non-Markovian. What 
we are going to develop is, on the contrary, a formula- 
tion valid in any case, also including our definition of 
non-Markovianity. 



A. Renormalized entropies 

Analogously to what has been done in the Markovian 
case, we start by analyzing the entropic ensemble behav- 
ior of the open quantum system. The equation for the 
time derivative of von Neumann entropy, Eq. (Fn), is of 
course the same, but now the time evolution of popu- 
lations is affected by memory effects. To see where ex- 
actly these effects come into play, let us closely analyze 
Eq. Q which follows from any TCL master equation, 
either Markovian or not. Such an expression for fib(t) 
depends on the transition rates Rb a (t) defined in |5| . No- 
tice however that each term in the sum on the right-hand 
side of ([5]) is proportional to a decay rate 7i(t) which, in 
the non-Markovian case, may cause the whole sum to 
temporarily become negative. To deal with this, let us 
rewrite the total transition rates as 



Rba(t) = Rbait) - R ba (*)) 



(19) 



where the Markovian contribution R\^,{t) and the non- 
Markovian contribution R^ a M (t) are defined by 

R-Z(t) = l[\R ba (t)\+R ba (t)], (20) 
RT(t) = l[\R ba (t)\-R ba (t)]. (21) 
With these definitions, Eq. ([7| can be rewritten as 

S(t) = -J2^(t)(Rl(t) - Inf£l (22) 

The occurrence of memory effects has the consequence of 
reducing the rate of change of entropy of the system dur- 
ing certain intervals of time. It is possible, as shown, to 
separate the non-Markovian contribution to the change 
of von Neumann entropy and to prove that such a con- 
tribution always counteracts the Markovian one. 

The next step is now to define entropy flux and entropy 
production for the analyzed process. In particular it is 
interesting to investigate the possibility of singling out 
non-Markovian contributions in these two quantities, as 
done for the von Neumann entropy itself. We then pro- 
ceed as before, writing 



S(t) = 6 e (t) + &i(t), 



(23) 



where formally both & e (t) and &i(t) are the entropy flux 
and production for the ensemble dynamics, and have the 



same expression as in the Markovian case, but they are 
fundamentally different because of the new structure of 
the transition rates Rf, a (t). Indeed, we have 



S e (t) 



S<(t) 



J2» 

a , b 



x In 



a (t)(Rli(t)-Rl M (t) 
Rfi(t)-RT(t) 



Ral(t)~RT(tr 



(24) 



a . b 



la(t)(RH(t) 



RT(t) 



x In 



„ (t) [Rm-RT(t) 



fib 



(t)(Rii(t) 



RTW 



(25) 



However these quantities, as already highlighted in [15] . 
involve logarithms of not necessarily positive terms and 
may become temporarily ill-defined (although both the 
von Neumann entropy and its time derivative are of 
course always mathematically well-defined) . Moreover, it 
is not possible to clearly isolate a Markovian and a non- 
Markovian contribution to these quantities, since the ar- 
guments of both logarithms involve Markovian and non- 
Markovian effects in a non-factorizable way. 

How is it then possible to speak about ensemble en- 
tropy production along non-Markovian dynamics? The 
appearance of negative rates is, by our definition, the 
characterization of non-Markovianity and, moreover, it 
is known that these rates cannot all be negative at the 
same time, which means there will always be at least one 
negative ratio involved in the definition of entropy flux 
and production. To overcome this problem, let us con- 
sider again Eq. (23). It simply amounts to writing the 



time derivative of von Neumann entropy as a sum of two 
contributions, each of which may involve the logarithm 
of a negative number. Notice however that, apart from 
the factor [i a j ^b (which can be written as an additional 
logarithmic term in the sum) , the arguments of each log- 
arithm in & e (t) and &i(t) are the same. If we replace 
all negative ratios R ba /Rab under the logarithms by their 
moduli \R ba / Rab\, the decomposition of S(t) into the two 
contributions still holds, since In \Rb a /Rab\ is added and 
subtracted in the sum. It is thus natural to define entropy 
flux and production for the non-Markovian ensemble dy- 
namics as 

S e (t) = - £>.(*) (j2g(t) - < M (t) 

a.b 



Rlt(t) 


-RT(t) 


Rai(t) 


-RT(t) 



x In 



a, 6 



x In 



(26) 



Ma(*) 




-RT(t) 


W>(*) 




-RT(t) 



(27) 
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These definitions coincide with the ones given in Eqs. ( 24 ) 



and ( 25 ) when the latter are real quantities, and ex- 



tend them to general non-Markovian ensemble dynam- 
ics. Equations (22 1 and (27) clearly show how memory 



effects manifest themselves in backfiows of information 
from the environment to the system as now Si(t) can 
become negative. 

There is however a second problem which is closely 
connected to the very definition of a fluctuation theo- 
rem: As in the Pauli master equation Q some rates are 
negative, it is not possible to give it a pure state single 
trajectory description. Nevertheless, it is possible to pro- 
ceed along a slightly different path which will indeed lead 
to a theorem for out-of-equilibrium fluctuations. To this 
end, consider Eq. ( 22 1 which is exact and directly stems 



from the TCL master equation we started from. It can 
be rewritten as follows, 



s(t) = -5>(*)(fl£(*) + j^RTit)) in 



= -^> (t)T bQ (t)ln 



a , 6 



Mb jt) 

Ma(*)' 



Mb ft) 
Ma ft) 

(28) 



where we have introduced positive renormalized transi- 
tion rates 



T ba (t)=K a (t) 



Mb ft) 



(29) 



Since positive transition rates characterize Markovian dy- 
namics we can exploit these renormalized rates to define 
effective Markovian-like flux and production for the open 
system entropy by means of 

Tba 



a,b 

^T(*) = Yl ^ aTba ln ■ 

^ Mb ^ab 

a, b 



(30) 



(31) 



These two quantities arc always well defined from a math- 
ematical point of view, and their sum just gives back 
Eq. (22 1. The microscopic motivation for these defini- 
tions will be given in Sec. |IVB| For now let us just look 
at Eqs. (30) and (31) as effective quantities, which on 



one hand solve the problem of negative arguments of the 
logarithms, and on the other hand reduce to Eqs. Q 
and (|9| in the limit of a Markovian dynamics. We can 
consider these quantities as the Markovian part of the 
expressions (26) and (p7|). The remaining part, which 



can not be effectively described as Markovian and which 
is thus irreducibly non-Markovian, is given by the differ- 
ence between the quantities (26) and (27), obtained as a 
direct extension of the Markovian ones, and the renor- 



malized quantities (30) and (31), 



S x (t) = S e (t) 



= ^ VaTba In 
a . b 



S r e (t) = Sl(t) 

Tba\Rab\ 



Si(t) 



\Rhn\T r . 



(32) 



We note that the quantity Sxft) is zero if all the renor- 
malized transition rates are equal to the original ones, 
i.e., if all transition rates Rba(t) are positive and there 
are no signatures of memory effects. 



B. Non-Markovian fluctuations 

In order to formulate non-Markovian fluctuations of 
physical quantities along single trajectories we first ob- 
serve that the master equation Q can be rewritten in 
terms of the renormalized transition rates J29l as 



Mb I 



ft) = ( T baftKft) - T ab {t)fl b (t) 



(33) 



This form of the master equation is strongly suggested 
by Eq. ( 28 ) for the time derivative of the von Neumann 



entropy and the corresponding decomposition into renor- 
malized entropy flux and production given by Eqs. (30) 
and pil. Note that Eq. p3l holds for both Marko- 



vian and non-Markovian processes, and that the renor- 
malized rates Tb a (t) are always positive by construction. 
As discussed in [24], the form (33) of the master equa- 
tion emerges if one interprets a negative rate i? oa for a 
transition from state a to state b as an effective posi- 
tive rate for the reversed transition which is given by 
Tba = r^lRabl according to Eq. Thus, we sug- 



gest employing the master equation ( 33 1 for the descrip- 



tion of fluctuations along single realizations of noncqui- 
librium processes. It should be noted however that in 
the non-Markovian case the transition rates Tb a (t) de- 
pend on the occupation probabilities and that, therefore, 
different trajectories are no longer independent which ex- 
presses the presence of memory effects [231 133 • 

Considering again a particular forward trajectory given 
by ( 10 ) we then find the corresponding probability 



N-l 



n*-n (t \ TT P -It' + 1 drJ2„T b (r) 

Pf - Mo i*o J 11 e J 

N-l 

xl[T a . +iaj (t j+1 )dt j+1 , (34) 
j=o 

simply by replacing the original transition rates by the 
renormalized ones. Correspondingly, the probability for 
the backward trajectory is given by 



N-l 



Pb - Majvl*JVj 11 e i 1 

3=0 
N-l 

x l[T ajai+1 (t j+1 )dt j+1 , (35) 

3=0 

and we obain for the logarithm of the ratio of forward 
and backward probability 



ln^=ln ^ (to) 



AT-l 



Pb 



T aj + 1 aj ftj + l) 



Ma w (*Ar) . Q T aja . +1 (tj + i) 



(36) 
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The first term on the right-hand side is again equal to 
the change of the von Neumann entropy As along the 
trajectory (see Eq. (15)). In analogy to the Markovian 
case, the second term represents the negative of the en- 
tropy flux integrated along the trajectory, i. e. we have 
(compare with Eq. (16)) 



N-l 



3=0 



T aj+iaj (tj+i) 

T aj a,j + 1 (tj + l) 



(37) 



Defining the renormalized single trajectory entropy pro- 
duction 



oy = As- As r e 



we immediately obtain from Eq. ( 36 ) 



P/(o>) 
Pb(-cy) 



(38) 



(39) 



Thus, we have found a fluctuation theorem for non- 
Markovian processes which is formally identical to the 
one obtained for Markovian dynamics, see Eq. ( 18 ). How- 



ever, in the non-Markovian case the fluctuation theorem 
holds for the renormalized entropy production ay which 
can be written as 



= As 



As 



As x = cr + As x , 



(40) 



where Asx is the single trajector y co ntribution to the 
time integral over Sx(t) (see Eq. (32)). In the Marko- 
vian case Asx vanishes and Eq. ( 39 ) reduces to Eq. ( 18 1 . 



This means that in general only a part of the entropy 
production, namely the one originating from the fluctu- 
ations described by the form ( 33 ) of the master equation 
can be described in terms of a fluctuation theorem. The 
fluctuation theorem for non-Markovian processes is thus 
fundamentally different from its Markovian counterpart, 
as it describes fluctuations of the entropy production of 
single processes, which are not the single trajectory con- 
tribution to the ensemble entropy production (27). 



To interpret the result (|39j) in physical terms we first 
note that within our approach fluctuations are described 
by the stochastic unraveling of the master equation (33 ). 



However, as has been emphasized already, due to the 
presence of memory effects single trajectories are not in- 
dependent of each other or, in other words, they are cor- 
related. Due to these correlations a single realization of 
a nonequilibrium non-Markovian process thus yields on 
average more information than just the one described by 
its associated entropy production. This is why the mea- 
sured renormalized entropy production, obtained in this 
work as ay in Eq. (40), is given by the usual single trajec- 



tory entropy production a plus an additional term which 
alters the usual Markovian form of the fluctuation theo- 



rem. This term, given by Asx in Eq. (40), represents the 



Correlations with 
other trajectories 



Single trajectory 
entropy 
production 




Net Markovian 
entropy production 



FIG. 1. (Color online) Schematic view of the contributions 
to the entropy producti on a long single nonequilibrium trajec- 
tories according to Eq. ( 40 1 . The net Markovian information 
flux oy, which obeys the fluctuation theorem (39 1, is equal to 



the sum of the single trajectory entropy production a and the 
non-Markovian contribution Asx due to correlations between 
trajectories. 



additional information extracted from the system origi- 
nating from the correlations between single trajectories. 
As we have demonstrated these contributions combine 



such that their sum obeys the fluctuation theorem ( 39 1 , 



completely analogous to the classical and the quantum 
Markovian one (see Fig.[7J). 



CONCLUSIONS 



The analysis of Sec. |IV| shows that, when taking into 
account memory effects in a quantum thermodynamics 
context, there is a close link between fluctuations in en- 
tropy production along single realizations of nonequilib- 
rium processes and the existence of an irreducibly non- 
Markovian entropic contribution in the ensemble dynam- 
ics. Any time the dynamics shows signatures of memory 
effects, a quantum fluctuation theorem has to take into 
account the full information contribution of the dynam- 
ics, which is no longer given only by the ensemble entropy 
production. More precisely, when a single trajectory is 
taken into account, memory effects produce an additional 
term to the information an external observer can reveal 
by measurements. Such an additional measured infor- 
mation contribution originates from the existence of cor- 
relations between single trajectories, which in turn is a 
consequence of memory effects. This can be clearly seen 
from Eq. (33), as the differential equations for the time 



evolution of populations are no longer linear due to the 
structure of renormalized rates, Eq. (29). Performing a 



measurement of entropy production along a single trajec- 
tory, then, means also extracting information about any 
other possible trajectory, and we have demonstrated that 
this extracted information leads to a net entropy produc- 
tion which effectively behaves as a Markovian one. It is 
this effective Markovian entropy production along sin- 
gle trajectories whose fluctuations can be described by 
means of a fluctuation theorem of the usual form. The 
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effective Markovian entropy production is obtained as the 
sum of the single trajectory contribution to ensemble en- 
tropy production and the information on trajectories cor- 
relations extracted by measurements. Quite remarkably, 
these two terms combine together in such a way that their 
sum behaves according to a very simple fluctuation law. 

Thus, due to non-Markovian features the measurement 
of fluctuating quantities affects the fluctuations them- 
selves. In a sense, the fluctuations we reveal in our ap- 
proach do not describe only single trajectory properties 
but supply information on the whole set of possible real- 



izations of a thermodynamic process and on their mutual 
dependence. The renormalization of rates performed in 



Eqs. (30) and (31 1 is, indeed, nothing but the mathemat- 



ical counterpart of the scheme depicted in Fig. [Tl We ex- 
pect, however, that the fluctuation theorem (39) might 



undergo more deep modifications if one considers other 
measurement schemes in order to characterize single tra- 
jectories, such as stochastic wave function unravellings 
expressing a continuous monitoring of the environment. 
This point could be an important and fruitful subject of 
future works in this field. 
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